Transport, optical properties and quantum ratchet effects for quantum dots and 

molecules coupled to Luttinger liquids. 
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We present non-perturbative solutions for multi-level quantum dot structures coupled to interact- 
ing one-dimensional electrodes out of equilibrium. At a special correlation strength the Hamiltonian 
can be mapped to the Kondo problem which possesses a solvable Toulouse point, where all con- 
ductance and noise properties can be calculated exactly. Special attention is paid to the fully 
asymmetric setup when each dot level is coupled to only one of the leads and the electron transport 
through the structure is accompanied by photon absorption (emission) . A relation between the opti- 
cal spectra and the energy dependent current noise power is established. Experimental implications 
of the results, specifically for the Fano factor, the ratchet current, and field emission via localised 
states, are discussed. In particular, we predict that the peak in the ratchet current as function of 
the irradiation frequency splits up in two due to correlation effects. 

PACS numbers: 73.63.-b, 71.10.Pm, 73.63.Kv 



I. INTRODUCTION 

Manufacturing of micro- and nano-electronic circuitry 
based on single molecules represents one prospective way 
to achieve further miniaturisation as well as efficiency 
improvement of electronic devices. First successful at- 
tempts of contacting single molecules have been reported 
in a number of recent experimental works [l|, y, Q . One 
possible mechanism for the electron transport through 
them is tunnelling on and off the molecular orbitals 
(MOs). The smaller the molecule the larger is the en- 
ergy distance between the MOs so that in some cases 
the transport occurs through only one electron level even 
at room temperatures. Hence, the adequate physical 
description of such systems coincides with that of the 
single-state quantum dot (QD): a fermionic level coupled 
to metallic electrodes (we shall also call them 'leads' or 
'contacts'). 

If one aims at small device dimensions one has to go 
for one dimensional (ID) electrodes. Promising can- 
didates for such wirings are the carbon nanotubes in 
their single- wall version (SWNTs) HH. However, truly 
one-dimensional electron systems cannot be described by 
the Fermi liquid (FL) model. No matter how weak the 
fermion interactions are, they cannot be taken into ac- 
count perturbatively. It is well known that in the low 
energy sector the interacting ID fermions constitute a 
universality class of Luttinger liquids (LL), which dis- 
play a completely different physics than the conventional 
FLs 0- As a consequence, the electronic degrees of free- 
dom in sufficiently thin SWNTs are also described by a 
generalisation of the LL model ^ _9J . This has also been 
confirmed experimentally [lol| . 

In the most of existing experiments the contacts be- 
tween the molecule and the leads are quite weak. The 
optimal operation of the future nano-electronic devices 
is, however, expected in setups with small contact re- 



sistances. That can be achieved only in systems where 
the current-mediating MOs undergo strong hybridisation 
with the valence bands of the electrodes. In the case 
when two MOs with different symmetries couple to dif- 
ferent leads electronic transport can only take place via 
emission or absorption of photons. Transport resulting 
from such 'optical' coupling can be distinguished from 
the background transport by its dependence on the laser 
irradiation frequency. In spite of vast amount of contri- 
butions dealing with QD's transport and optical proper- 
ties (see, for example, |Tll |. and references therein), sys- 
tems with good contacts did not receive much attention. 
Moreover, to the best of our knowledge, the issue of in- 
teracting electrodes has been discussed only in resonant 
tunnelling context [H [H El El E 113 . In this paper 
we begin to close these gaps and present non-perturbative 
solutions for multi-level dots contacted by interacting ID 
electrodes. We wish to clarify here that these type of cor- 
relations are different from 'on-dot' couplings (Coulomb 
or Hubbard terms). The latter coupling gives rise to the 
Kondo type phenomena E3l- this paper we neglect 
the on-dot couplings and effectively deal with spinless 
electrons. This restricts the validity of our results to the 
temperatures above the Kondo temperature Tr- or to the 
case when a polarising magnetic field is present. Such an 
approach appears to be justified for SWNTs, which dis- 
play strong LL correlations in a wide temperature range 
from about 5 K to about 100 K E3|j typical T/^-'s being 
much smaller. The interplay of on-dot and 'in-lead' cor- 
relation is of theoretical interest but remains outside the 
scope of this paper. 

The outline of the paper is as follows. 

In Section ^ we start with a review of the simplest 
realisation: a single-state dot coupled to ID electrodes. 
These results have been recently announced in our Letter 
[13. Contrary to [13 we work in the Green's function 
formalism, which is more suited to access the noise prop- 
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erties. Since the transport in such a setup is completely 
understood both in the resonant case as well as in the 
off resonant case, our primary goal in the Section ^1 is 
the study of the Fano factor, which is a ratio of the zero 
frequency noise and the transport current, as a function 
of bias voltage. For the sake of completeness we derive 
all equations needed to access the full energy dependence 
of the noise power spectrum. 

In Section lTTll we go over to a two-state QD with both 
levels coupled to ID leads. Since the non-interacting situ- 
ation is more or less trivial we restrict our considerations 
only to the interacting system. It turns out, that, as 
in the case of the single-state QD, the Hamiltonian can 
again be brought to a quadratic form in terms of new 
fermions for a special interaction strength In this 

representation the non-linear transport as well as zero 
frequency noise properties can be expressed via the trans- 
mission coefficient even in case of an additional tunnelling 
term between the dot levels. Our non-perturbative ap- 
proach allows then to study all resonant tunnelling effects 
(known for single-state setups j3) ™ this situation. 

In Section IIVI we investigate transport in a similar 
structure where every level is only coupled to one of 
the electrodes and where the dominant transport mecha- 
nism is the photon-assisted tunnelling between the levels. 
That situation corresponds to coupling of two MOs of dif- 
ferent symmetries discussed above. It turns out that a 
finite current can flow even without any applied voltage. 
That makes such a system one of the simplest realisation 
of the so-called 'quantum ratchet' effect 19, 20, 21] . We 
concentrate on the analysis of the 'ratchet' current as a 
function of the radiation frequency. Some general results, 
including an important relation between the current noise 
power and the absorption and emission spectra are dis- 
cussed in Section FlV Al A treatment of non-interacting 
systems follows in Section IIVBI Transport through a 
dot coupled to LL electrodes is then analysed in Section 

A short summary of results (Section 01 concludes 
the paper. We stress again that though we model 
the molecule-electrodes coupling by means of tunnelling 
Hamiltonians, all our results are non-perturbative in tun- 
nelling amplitudes contrary to the bulk of existing studies 

m 

II. TRANSPORT THROUGH A SINGLE-STATE 
QUANTUM DOT 

First we briefly review the method of and then 
discuss the noise properties. 

A. Scattering states solution and duality 

We model the system by the following Hamiltonian (we 
ignore the spin throughout the paper): 

H ^ HK + Ht+Hc, (1) 



where Hk is the kinetic part, 

HK = AdU+ J2 Horn, 

i=R,L 

describing the electrons in the leads iJo[^i], and the reso- 
nant level A, the corresponding electron operators being 
The dot can be populated from either of the two 
leads {i = R, L) via electron tunnelling with amplitudes 

i/t -^7,[dt^,(0)-f h.c.]. 

i 

In iQ, He describes the electrostatic Coulomb interac- 
tion between the leads and the dot, 

i 

This interaction is a new ingredient we have introduced, 
absent in the related studies and 0|. It does not, 
however, affect the universality as we shall show later. 

The contacting electrodes are supposed to be one- 
dimensional half- infinite electron systems. We model 
them by chiral fermions living in an infinite system: the 
negative half-axis then describes the particles moving to- 
wards the boundary, while the positive half-axis carries 
electrons moving away from the end of the system. In 
the bosonic representation i7o['0i] are diagonal even in 
presence of interactions (for a recent review see e.g. 1^; 
we set the renormalised Fermi velocity v = vp / g = 1, the 
bare velocity being vp)- 

Horn = (4^)-^ J dx [d^M^)]''. (2) 

Here the phase fields 4>i{x) describe the slow varying spa- 
tial component of the electron density (plasmons), 

ipl{x)il;i{x) = dx<t>i{x)/2Ti^ . 
The electron field operator at the boundary is given by 




V>,(0) = e'^-(°)/v^/V2^, (3) 

where oq is the lattice constant of the underlying lat- 
tice model. Here g is the conventional LL parameter 
(coupling constant) connected to the bare interaction 
strength J/ via g = (1 4- U /ttvf)~'^^'^ Hill- In the chiral 
formulation the bias voltage amounts to a difference in 
the densities of the incoming particles in both channels 
far away from the constriction p^ . The current is then 
proportional to the difference between the densities of 
incoming and outgoing particles within each channel. 

To the best of our knowledge, Hamiltonian |^ cannot 
be solved exactly even in the g = 1 case as long as Ac 
remains finite. However, after a transformation of rf^ and 
d operators to the spin representation of the form 

= (d^ + rf)/2, 
Sy = -iid'^ -d)/2, , 
= dU-l/2, 
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one immediately observes that the Ac term is analogous 
to the S^-spin density coupling in the Kondo problem. 
The latter is known to be explicitly solvable at a particu- 
lar value of the longitudinal coupling: the Toulouse limit 
(see e.g. 0). Let us perform a similar calculation. As a 
first step we introduce new symmetric and antisymmetric 
fields 



±1 



(4) 



which still fulfill the bosonic commutation relations. 
Then we apply the transformation H' = WHU with 

121 



which changes the kinetic and the Coulomb coupling 
parts of the full Hamiltonian to [we drop a constant con- 
tribution proportional to SzS{x) that does not affect the 
transport] 



and the tunnelling part (terms containing 7.;) to 



where S± — Sx iSy — d\d. At the point g = \/2 one 
can re-fermionise the problem by defining new operators 



(5) 



which fulfill standard fermionic commutation rela- 
tions. With the help of the particle density operator 
V'-LV'i — dx4>±/'2TT we can immediately write down the 
refermionised Hamiltonian, 

H = Ho[i>±] + iXc~2n)2S,^l^P+ + AS, 

+ 5'+(7lV''- + + ilLtp- + 'yRil;^)S^ . (6) 

In the case of the symmetric coupling 7^ = this 
Hamiltonian is similar to that of the two-channel Kondo 
problem and, at the Toulouse point Ac — 27r, can be 
solved exactly (out of equilibrium) using the method of 
Ref . [231 . The novel ingredient in the following analysis is 
the extension to the asymmetric case. To take advantage 
of the Toulouse point we set the Coulomb coupling am- 
plitude to 27r in what follows. This not only removes the 
four fermion interaction but decouples the '±' channels 
making the '+' channel free (i.e. decoupled from the dot 
variables). 

At the Toulouse point our Hamiltonian describes free 
fermions which scatter at the origin. These new fermions 
are related to the physical electrons in a highly non-local 
way. As the relations between the particle densities are 
still linear, in order to access the transport properties it is 
sufficient to calculate the energy dependent transmission 



coefficient 1 — T{uj) of the new fermions [T{u!) being the 
transmission coefficient of the physical ones] . The non- 
linear I — V characteristics is then given by 

I{V)^GoJ du;T{u;)[nFiio-V)-nFiuj)], (7) 

where np denotes the Fermi distribution function and 
Go = e'^/h is the conductance quantum. The easiest way 
to identify T{uj) is the equation of motion method. We 
calculated T{u!) in Ref. [la, it is given by (we measure 
all energies in units of F = 7^ -(- 7I,): 



(8) 



{E^ + PDiE^ + 131) + 2j^{E^ + /3_/3+) + 7^ 



where 



E = A^-u;\ 
/3± = [(l-2a)A±c^]/2, 



7 = u!\/a{l - a) , 

and a — 7£/(7|, -I- 7^;) is the asymmetry parameter. 

Using expressions JTJ and ||SJ) one can access all con- 
ductance properties of the system. We shall not discuss 
them again (see Ref.il^]) but concentrate instead on the 
duality property. In the simplest case of the symmetric 
model on-resonance {a — 1/2, A = 0) we obtain 



T{lu) = (l+cj2)- 



(9) 



As a consequence, the temperature dependent differen- 
tial conductance at zero bias (which is the most relevant 
quantity from the experimental point of view) amounts 
to 



GA=o(r)/Go 



1 



2ttT 



1 



1 



2ttT 



(10) 



where ip denotes the ^-function. Comparing this result 
with the conductance G1/2 through a single scatterer in 
an LL with g ~ 1/2, given in Ref. [3], we find that 
(similar relations are valid for the non-linear / — V^'s) 



Ga=o(T)/Go = 1-Gi/2(T)/G, 



(11) 



where T in Gi/2{T) is measured in units of the backscat- 
tering strength. According to the duality hypothesis the 
strong coupling fixed point at g corresponds to the weak 
coupling one at 1/g and vice versa, leading to the rela- 
tionship of the conductances of the form [l3. 125I |2^ , 



G2(T)/Go = l-Gi/2(T)/G 



(12) 



Therefore, Eas. ((TT)) and l(T^ suggest that tunnelhng be- 
tween two LLs with g — 1/2 via a resonant level is equiv- 
alent to direct tunnelling between two LLs with g = 2. 

In fact, we can demonstrate this equivalence on the 
Hamiltonian level. We start with ® and in order to 



4 



simplify things introduce new real (Majorana) fermions 
a, b and C, rj according to 

d^{a + ib)/V2 , V- = (C + iv)/V2 . (13) 

In this language the Hamiltonian acquires the form 

H = Hole, V, a, b] + i7_a77(0) + ij+bm , (14) 

where j± — 7^ ±7/j and the unperturbed part is defined 
by 



HQ[(,ri,a,b] = iAab + i / dx 

+ vcixMx) 



For future reference, we give the current operator in this 
representation: 

J = -^7-aC(0)-^7+^(0). 

On the other hand, the Hamiltonian for the direct tun- 
nelling between two LLs in terms of the physical fermions 
is 

H= J2 Horn + 7 (^iv-fl + v-Ul) , 

i=R,L 

containing the free part Hq, which again describes two 
half-infinite LLs, and tunnelling between them with the 
amplitude 7. We bosonize the above Hamiltonian as in 
the previous Section using rules ^ and (|3J), and obtain 
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27rao 



Introducing new fermions © at the point g = 2 (we drop 
the channel again since it is free) 



H = Ho[iJ-] 



Next we take advantage of a trick from Ref.[23 and per- 
form a substitution 5* = (d — d^)ip^, where d is some 
local fermionic operator not related in any way to the 
dot operator of the previous Section. Obviously, such 
transformation does not change either the commutation 
relations or the normalisation of the operators. Then 
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\/27rao 



(c-c^)(1'^ +*), 



The last step is obvious: one introduces the Majorana 
components according to (|13|) . This results in 



H 



TTflo 



which is precisely the Hamiltonian H14(l of a resonant 
setup (7- = 0, A = and, of course ^ = 0) up to the 
redefinition 7+ = jy^2/nao. We have checked that the 
bias voltage and the current operator of the g — 2 LL 
problem and the resonant tunnelling system transform 
correctly. 



B. Green's functions solution and noise properties 

Although the transport properties can easily be ac- 
cessed by means of the scattering formalism as shown in 
Ref . [1^ ■ it is not immediately clear (see below though) 
how the information about the finite frequency noise can 
be extracted from the transmission coefficient. A more 
appropriate method to calculate the fluctuations is the 
Green's functions (GFs) method in its non-equilibrium 
(Keldysh) formulation. For further reference we now de- 
fine all possible non-equilibrium GFs. Let /i, v stand for 
either of the electrode Majoranas C or 77 (taken at x = 0) 
and /, h stand for either the dot level operators, a or b. 
Then define 



D}\{t-t') = 
G%{t-t') = 



-i{Tcf{t)Ht')) , 
-i{Tcpi{t)v{t')) , 
-i{Tcii{t)f{t')) , 
-i{Tcf{tW)) 



where Tc is the time ordering operator along the Keldysh 
contour C, which consists of the forward C_ and the 
backward C+ paths. The times t' belong to the paths 
Ci_j, respectively. Sometimes we shall omit the Keldysh 
indices i,j in what follows adopting matrix notation for 
the Keldysh GFs. In the above definitions we assumed 
the system to be in a steady state so that all GFs are 
translationally invariant in the time domain and there- 
fore depend only on the time differences. There is an 
obvious relation G^^j{t — t') — ~Gj^{t' — t). However, 
working with two definitions possesses advantages as we 
shall see later. 

It is, in fact, not difhcult to calculate the zero order 
GFs, when 7^ = 0. For the electrode Majoranas we ob- 
tain 



g[?(c.) = G(°)(c.) = 1 



H{lo) + 1 
H{u:) 



(15) 



where H(lo) contains the information about the Fermi 
distribution functions of the original electrons in the 
leads, H{ijj) — np{uo + V) — np{—ui + V). Obviously, the 

cross correlations G^°^ exist only as long as the applied 
voltage is finite. 



Gf^{u^) = -G':'>{u) = ^F{u^) 



(0), 



'C,v 



1 1 
1 1 



(16) 



This fact is reflected by the function F{u!) = np(u) + 
V) — np{u! — V) vanishing as V 0. As a consequence 
of the special forms of (|15|l and 1)16(1 the retarded and 
the advanced components are fairly simple: G^°2^'"* — 



while G 



{0)R,A 



=pj/2. The a — b subsystem being in 



equilibrium makes the calculation of the corresponding 
GFs even simpler. The result is 



D 



(0)^-(++) 

if 



0=± 



5 



Dfj-+^+~^ = ±i7r5{u + A), 

%=± 

Using these zero order GFs as a starting point we can 
calculate any correlation function exactly because of the 
quadratic form of the Hamiltonian. Our goal is to calcu- 
late the average of the current operator I{V) = (J) and 
the power spectrum of current fluctuations (we also shall 
call this quantity noise spectrum), which is defined as 

P{n) = I dte^"' [(J(t)J(O)) - (J(0)>2] . (17) 

The averages are calculated using the S matrix for the 



coupling of dot and electrode Majoranas, (...) = {. . . S)o, 
which is 

S = Tcexp ^- J dT-f_a{T)ri{T) + 7+6(t)C(t)^ . 

Expanding in powers of one can derive for the current 
the analog of the Meir-Wingreen formula , 

W=8^ / dc.F{u;)[^lDt,iu;)^^lDtic.)]. (18) 

We choose to split the noise spectrum into two contribu- 
tions, /'(ri) — P\\{^) +P±{^1), where both quantities can 
be expressed in terms of off-diagonal (containing different 
fermion species) GFs, 



Plim = - I dco^i G+-(^)G+-(f7-c.) 



p±m = 



7-7+ 



dcu 



- D+-icj)G+-in-Lu) 

+ G+- {uj)G+- (Q-uj)- 2D+- {u)G+- {n - lo) 



(19) 
(20) 
(21) 



By means of the S matrix expansion one can reduce 
some of the off-diagonal GFs to the diagonal ones. In 
particular, for the functions entering the parallel part 
of the noise spectrum we obtain (we omit the energy 
variable w), 

Gt = ^\FDg - ^l-{D^T-G^, + dTg^,;) (22) 

Similar expressions can be derived for G^j^ , G'^^ and 
G't~ . GFs entering the Pj^ part of the noise have some- 



what different structure, 

where the remaining GFs, G'^^ , and G^j~ are sim- 
ilar. The remaining off-diagonal GFs cannot be reduced 
to the diagonal ones. They should rather be found as 
solutions of a chain of Dyson equations. For the a — 
subsystem one obtains the following system of equations 
(to simplify notation we ignore here the Keldysh indices) : 



Da, - ^i".^+7i^i°;Gg)i?., + 7^i?i"jG,,i?i°)+7+7-(i?i°)Gf>.^ (23) 

D,, = Di"J + .iD^^Gf^D,, + 7^i?I°^G(o)Z)., + 7-7+(^I°)g(? A. + (^4) 
In the same way one can derive the corresponding equations for the electrode Majorana GFs: 

Gc, = G[; + jlGf^^Dil^G,, + 7^G(>i"jG,, - j^j^{Gf^^ dI'Jg,, G^^^dH^G,,) , (25) 

G,, = G(°) + 7iG(?i?i°)Gc, + llG(^^D(:jGr,r, + 7+7- (Gi?i?i°jG,, + G^^^D^^Jg,,). (26) 



In the simplest symmetric case 7_ = we obtain from Ea. (|24(l for the advanced dot level GF 
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Plugging this result into the expression for the current 
(|18|1 results in Eq.Q with all the energy variables mea- 
sured in units of F = 7^ /4 we identify 



- LO 



(27) 



as precisely the transmission coefficient ((SI at a = 1/2 
found previously by means of the equations of motion 
method. 

As for the noise spectrum, only the parallel component 
survives. Ea. H24l) has the following solution: 



D 



,(0)- 

'bb 



,(0)i?|2g(0) 



bb 



p^(0)«^(0)«|2 



while its electrode counterpart, 



Q + - ^ q{0) + - 
?7?7 rjrj 



is related to the cross-correlation G^^. 
solution to one of the equations in (|25|l . 



The latter is a 



Cv 



which is zero. 



Therefore we have G^^ 



G 



(0)+- 



The 



rest of the needed GFs can be read off Ea. l22|l . G'l^ — 
ij+FDg/2, and G+- = -n+FD^j2. Collecting all 
terms in (|19|l we obtain the following result: 



Psymi^) = eGo duj \ - 



F{Lu)F{n~uj)Lu{n-uj) 



I 



[g(cc;)-l][g(»-Lj)-l]c 



(28) 



The same formula has been obtained by Schiller and 
Hershfield (SH) in Ref.l?! in the context of the non- 
equilibrium Kondo problem, where the magnetic field 
strength plays the role of our level energy A. 

In the following we do not repeat the results for the 
noise spectrum which are already contained in Ref . |24| 
but rather concentrate on the aspects which are specific 
to resonant tunnelling, in particular on the asymmetric 
case and on the calculation of the Fano factor not covered 
by SH. 

In the limit of zero frequency we find the noise spec- 
trum to be given by the formula 



%,„(0) = eGo / <hjT(u:){\^T{i^)\ 



X \np(io — V) — np(uj)\ . 



(29) 



identical to the one derived for non-interacting electrons 
[2 9*1. This is somewhat surprising. The reason must be 
that we map the original Hamiltonian onto a free one, 
where the current carrying excitations are again of the 
fermionic nature. That is why in order to access P(0) 
even in the asymmetric case we do not have to solve the 
above Dyson equations but can simply use formula (jS)) 
for the transmission coefficient. 

Let us pause here to mention that, in non- interacting 
resonant tunnelling systems at zero temperature, the 
Fano factor vy = P{Q)/eI{V) is suppressed in com- 
parison to the Schottky value z/ = 1 by the factor 
{Tl+T\)/{VL+TRf at high voltages F > F and by the 
factor {T L-T / {T L+T in the opposite limit F < F 



[33. As previously Yl{r) denote the dimensionless con- 
ductances of the left (right) contact. The suppression is 
maximal, 1^00 = 1/2 and vq = 0, in the symmetric case 
Tr = Tl. According to Ref.lU, as soon as we deal with 
an LL system, the Fano factor is expected to keep its 
maximal value no matter how strong the asymmetry is. 

At zero temperature and on-resonance we find using 
^ that the on-resonance shot noise is (we also recover 
the correct pre- factors): 



^sj/m(0) = [tan 



V - V/il + V"^ 



Taking into account the formula for current, I{V) = 
Gotan"^ y, one can read off the Fano factor, which has 
the following limiting forms: 



(30) 



In the general asymmetric situation the transmission 
coefficient takes the form 



T{u;) 



[cj2 + 1/4 + 0(1 - a)\ - a(l - a) 
The evaluation of the Fano factor now yields 



2a^ - 2a + 1 



p2 , p2 



(31) 
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FIG. 1; The Fano factor as a function of the bias voltage for 
different asymmetry values (a — 0.5,0.4,0.3,0.2,0.f from the 
bottom curve upwards). 

At high voltages, V ^ T, we recover the noise suppres- 
sion of the non-interacting case. On the contrary, our 
value for Ml is in apparent contradiction to the results 
of Ref. [31|. The reason for that discrepancy is quite 
simple. Ref. [3lj | assumes the sequential tunnelling pro- 
cess to be the dominant transport mechanism. As was 
pointed out in 0, this is indeed the case for not too 
low temperatures for arbitrary g and even at T = (that 
is in our situation) as long as the interactions are strong 
enough, for < g < 1/2. These conditions are obviously 
not compatible with our assumptions. 

It is, in fact, not difficult to access the full crossover 
behaviour of the Fano factor, see FigQ The most strik- 
ing feature of the full plot is the presence of the local 
minimum at V* as long as the system is kept either sym- 
metric or out of resonance. V* can be shown to be the 
solution of = T{V*), so that it gives precisely the 
point at which the transmission coefficient crosses uy- 
Similar local minima have also been found in [Slj. 

In the case of a symmetric system off-resonance, where 
the transmission coefficient is given by Ea. H27|l . the 
emerging picture is fully consistent with H31|l . In the 
small bias limit the Fano factor approaches unity, vq — 1, 
whereas in the limit of high voltages we again recover the 
universal non- interacting noise suppression as v^o = 1/2. 
The vv behaviour is qualitatively the same as in the 
asymmetric case, see Fig. ^ including the minimum at 
intermediate voltage. The asymptotic value at Uao = 1/2 
is now universal for all curves. 

We performed a study of the general case a ^ 1/2 and 
A ^ as well. The limiting behaviour of vy turns out 
to be determined solely by the asymmetry parameter a 
and is completely independent of the detuning A. For 
A ^ 1 the position V* of the intermediate minimum is 
asymmetry independent and coincides with A. 

Contrary to the zero frequency noise the evaluation of 
P(r^) spectrum at finite O and, in the general case of an 
asymmetric system off-resonance, requires knowledge of 
the full transmission amplitude matrix, as one already 



can see from ^ [a^. The latter is not at all the same 
in both fermionic representations of the problem, so that 
we have to solve the full set of Dyson equations 
The full solution is rather lengthy. It does not appear to 
reveal qualitativel y n ew features as compared to what 
is already known [2J| as the effects of asymmetry and 
finite A do not compete but rather enhance each other. 
We shall therefore conclude the discussion of the noise 
spectrum at this point. 



III. THE TWO-STATE QUANTUM DOT 

In the spirit of the previous Section we model the dou- 
ble QD by two fermionic levels with energies Ai^2, which 
are coupled to LL leads, see Fig. |21 As long as the levels 
do not interact with each other in any other way then by 
tunnelling, the whole treatment including the solution of 
the equations of motion can be performed for an arbi- 
trary number of levels. Throughout this Section we are 
not interested in noise power spectra so that we concen- 
trate only on the conductance properties of the system, 
which are most easily accessed by means of the equation 
of motion method. The Hamiltonian of the system is still 
assumed to be of the form |^ with following changes: 

(i) The kinetic part describes two (or more) levels in- 
stead of only one. 



E 

i=l,2 



A,- did. 



E 



(ii) The tunnelling amplitudes are the same for both 
levels. It is, in fact, not difficult to solve the problem with 
arbitrary amplitudes. This, however, does not induce 
new physics, so we restrict our solution to this special 
case: 

^* = E E 7,[4^.(0) + h.c.]. 

i j=R,L 

(iii) The strength of the electrostatic Coulomb inter- 
action is also assumed to be the same for both levels, 



i j 

(iv) There is an additional term in the Hamiltonian 
that is responsible for the tunnelling processes between 
the dot levels Q: 

Hw = W{d\d2 + h.c.) . 

To describe the electronic degrees of freedom in the 
electrodes we use the same formalism as in the previous 
Section, see Eqs. As in the case of a single level 

we can introduce symmetric and anti-symmetric compo- 
nents Q and spin representations for the level operators 
(they acquire an index z = 1, 2) and apply a slightly dif- 
ferent transformation to the overall H, defined by 



f/ = exp z^5,^0+/V2^ 



(32) 
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FIG. 2: Quantum dot with two levels. 

This transformation changes the kinetic and the Coulomb 
coupling parts of the full Hamiltonian to (we again drop 
a constant contribution) 

H'k + H'c = Hk + {Xc/tt^ - y^) ^ 5,^9,0+ (0) , 

i 

and the tunnelling part (terms containing 7^) to 



where — Sf ± iS^ — dl,di. The intra-dot tunnelling 
term Hw is invariant under this transformation. The 
refermionisation can again be performed using the def- 
initions (O and the resulting Hamiltonian differs from 
that in ((BJ only by the sums over both spins, 

i 

+ SfiiLi^- + iR^l) + (7L^i + 7fl^-)^r -(33) 

In what follows we concentrate on the Toulouse point 
when Ac = 27r, where the '±' channels decouple 
and when the Hamiltonian acquires a very convenient 
quadratic form. 

In order to calculate the non-linear I{V) we employ the 
method of Ref . [T^ . which results in Eq.jTI) with some 
modified r(w). As in the case of the single-level dot, 
the easiest way to find the transmission coefficient is the 
equations of motion method. Since we have two types of 
operators: for the electrons of the '— ' channel and for the 
dot levels (we go back to the original d\,di operators), 
we need two types of equations of motion, 

idtip-ix) = -idxi/j-jx) +y^^S{x){'-fLdi - ^Rdl) , 

i 

idtd, = /\dr+Wd^,+-iL4'-{^)+lR^lmM) 

Integrating the first one around x — Q with respect to x 
from — e to e and then sending e to zero we obtain 



where 0^ denotes positive (negative) infinitesimal. We 
define new operators 

i 

L = 2{idt + VK) - ^ A, , 

i 

Z± = ^■^ + iY,^^^t±{AoAl-W^) . (36) 

i 

By acting with \Y\'^ on both sides of Ea. (|35|l and using 
the last two equations (|34|l we can eliminate the dot op- 
erators. We obtain as a result 

|ynV'(o+)-v(o-)] = hlz+L + j%z^L*)m 

+-fLjR{Z+L + Z^L*)i;^Q). (37) 

Now we can insert into this relation the momentum de- 
composition of the field operator ijj^ 

*-(M) = /|.-<'-"{-t">%° ■ ™ 

Because the dispersion relation is linear, lo = vk = k, 
we can use as the momentum variable as well as the 
energy variable. Inserting Ea. (|38|l into Ea. (|37|l and us- 
ing ip-{0) — [tp-{0~^) + ip-{0~)]/2 results in a following 
equation, which form (up to a redefinition of constant 
factors) is independent of the number of levels, 

E{b^ - a^) = -iP+{a^ + b^) + ij{al^ + blj . (39) 
We introduced the following objects. 



E = 



7 = 



n 



p— ± L i 



Y[{Lu+pAi)-W 

i 

v/a(l -a)[^ A] - 2uj 

j 

2W A, - 



/3± = (I/2-a)[^A, n(A^-^')-2W 



H^2^aJ ±7/2Va(l-a) 



*[^_(0+) - V'-(O-)] = V7Lrf. - 7fl4 



(35) 



(40) 



Formally Eq. (|39() has exactly the same form as in the 
single-level case. Therefore the resulting transmission co- 
efficient is still given by formula ijHl with modified con- 
stants contained in (|40|l . 

The temperature behaviour of G at the maxima does 
not differ considerably from that for one single level, 
which has already been studied in Ref. . The conduc- 
tance G{T) in the valley between the peaks turns out to 
be very well described by a superposition of two peaks of 
the single level problem. Moreover, the presence of the 
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FIG. 3: Linear differential conductance of a symmetric setup 
as a function of the gate voltage Vg at the temperature T — 
O.Oir at different values of the intra-dot tunnelling amplitude. 
The bare dot level energies are Ao = — F and Ai = 0. 
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FIG. 4: The same plot as in Fig|21 but for an asymmetric 
system with a = 0.25. 

W tunnelling process does not affect the conductance ei- 
ther, only leading to some suppression of current at low 
temperatures T ^ Tj^ l. 

The issue of the peak shape is far more interesting, 
see Figs. I3I4I In the asymmetric system the upper peak 
tends to sharpen and to increase in height with growing 
tunnelling amplitude whereas the lower peak suffers the 
opposite fate, see Fig. 01 This can be understood in 
the picture where the Hamiltonian is diagonaliscd with 
respect to W tunnelling. The corresponding transforma- 
tion is given by a rotation in the two-dimensional space 
of operators di,2 to the new ones di_2 via the rotation 
matrix R{a), where the angle a is given by 



tan 2a 



2W 



A2- Ai 

Then the energies of the two new levels are 



Ai,2 = 



Ai A2 



± 



'Ai- A2 



The couplings to individual levels are subject to change 
as a consequence of the operator sum transformation 

di +d2 = \/2cos(a + 7r/4)fii + \/2cos(q; - n/4)d2 , 

which makes the coupling of the upper level decrease 
with growing tunnelling amplitude W when a tends to 
7r/4. The coupling of the lower level, on the contrary, 
increases. Notice that such renormalisation occurs even 
in the non- interacting systems. The physical reason is 
that the conductance via tunnelling through the lower 
level is enhanced because of the additional depopulation 
process which transfers electrons to the upper level. The 
non-trivial interaction effect in our setup is the increas- 
ing height of the upper level. The reason is that due to 
smaller Tu l the upper level is at effectively higher tem- 
perature, which, in turn, means that the conductance is 
higher in the asymmetric case. For the same reason the 
height of the lower peak diminishes. 

The coupling symmetry is not affected by the diagonal- 
isation transformation, so that the resonant conductance 
is increasing monotonically all the way to zero tempera- 
ture. That is why in this case the upper peak amplitude 
is lower, see Fig. 13 



IV. QUANTUM DOT INTERACTING WITH 
PHOTONS: A QUANTUM RATCHET SETUP 

A. Preliminary considerations 

Now we slightly change the setup. The kinetic and 
Coulomb coupling terms Hk+Hc remain the same while 
the coupling to the leads becomes completely asymmet- 
ric: the level Ai is coupled only to the left lead and A2 
only to the right lead (see FiglS}, 



(41) 



We assume that the localised levels possess different sym- 
metries so that a direct tunnelling between them is for- 
bidden (nevertheless, the Coulomb coupling, being free 
of selection rules, is still symmetric) while hopping with 
simultaneous emission or absorption of a photon with en- 
ergy is possible. Then the coupling of the levels can 
be written as 



H 



w 



W 



e'^'d\d2 + e 



dkdi 



where W denotes the coupling amplitude. We dropped 
the photon operators since we assume that the dot (or 
molecule) is subject to intensive laser radiation (so that 
there is always a phonon which can be absorbed). 

Although the full Hamiltonian is now explicitly time 
dependent it still can be reduced to a system in a steady 
state via the gauge transformation 



(42) 
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FIG. 5: The quantum ratchet setup. 

Thereby the energy levels of the dot states as well as the 
chemical potentials in the leads are shifted to 

Aia Ai^2Tf^/2, 

Thus, the coupling to the radiation effectively results in 
a finite bias voltage. As a consequence, a finite cur- 
rent can flow without any real applied voltage at all. 
This is the quantum ratchet effect which has recently 
been intensively discussed in Ref . 0, [23, |^ in various 
non-interacting setups with mostly weak couplings to the 
leads. 

The physical explanation of the effect in the original 
picture prior to the gauge transformation is quite simple. 
We start with a system where the two dot levels possess 
two different energies, A2 < Ai as in FigO The popula- 
tion probability of the lower of level (A2) is higher than 
that of its counterpart Ai, so that the photon absorp- 
tion is the dominant process transferring the electrons to 
Ai. They can relax either to the left lead or back to A2. 
However, if the hybridisation of the dot levels is larger 
than the electromagnetic coupling, the dominant relax- 
ation process is tunnelling into the left lead. This leads 
to a non-zero net current through the system. 

The quantities we are interested in are again the full 
current I{V) through the system, the current noise power 
P{uj) and the light absorption (emission) A{lu)[E{lu)] 
spectra. The first quantity can be defined e.g. via the ex- 
pectation value of the particle flow between the left lead 
and Ai, 

i=-ijL[d\MO)-^liO)di]. (43) 

In the language of non-equilibrium Keldysh diag ram 
technique it is given by the following expression |33l l3^ 



2tt 



duj 



G 



(0) 



-(^)i?r+(^) 



G 



(0) 



-(^)I?+-(^ 



(44) 



where — -f (H — ) indices stand for lesser and greater 
Keldysh GFs. For the left (right) lead electrons they 
are denoted by G^^(a;) while for the dot electrons by 

D\^2{^)- The additional superscript (0) distinguishes the 
GFs in the absence of tunnelling couplings. Alternatively 
one can define the current as the expectation value of the 



tunnelling operator between the dot states, 
/' = -iW{d\d2 - dldi) . 



(45) 



The third possible method to describe the transport is 
the transmission coefficient formalism. Using definition 
(I45|l we derive an expression for the noise power spectrum 
of the form: 



dte^' 



{d\{t)d2{t)dl{0)di{0)) 



(4 (i)di(t)4 (0)^2(0)) - 27TS{uj)i 



The light absorption rate can be evaluated via the 
Golden Rule and is given by 

A{lj) = 27rVK2|(/|4d2|0)p5(c^- Ai + A2), 

where the vacuum |0) is assumed to be the state with the 
level Ai empty and A2 full. The final state |/) is just 
the opposite. It turns out that the Fourier transform of 
the function 



(46) 



S^it) ^ -i(dl{t)di{t)d{{0)d2{0)) 
is directly related to the absorption rate, 
A{lu) ^tW^S^iiu). 



Another interesting aspect of Eq. (|46() is the fact that it 
can be written down in the form 



Ztt 



deD:^+{e)D+-{e + Lu). 



(47) 



This formula is exact as long as the GFs are calculated 
exactly and the level operators participate only in two- 
particle interaction terms (such as Hw or Ht). 

Similarly one can show that the emission rate E{uj) is 
proportional to a related function 



where 



S^{t) = -i{d{{t)d2{t)dl{0)di{Q)) 



(48) 



(49) 



the Fourier transform having the same form as (|47|) up to 
the exchange 1 2. Thus we establish a very convenient 
expression relating the total optical spectrum with the 
noise power of our system: 



P{lu) ^ A{lu) + E{lu) - 2tt5{uj)P 



(50) 



This relation has far reaching consequences for the field 
emission (FE) physics In the case when the chemical 
potential of the right lead is sent to fin = — 00 the dou- 
ble dot setup describes tunnelling into vacuum through 
a sequence of localised states, which is nothing else than 
FE |37l |. Experimentally, it turns out that under certain 
conditions the FE from carbon nanotubes is accompanied 
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by luminescence phenomena 1381 . One possible explana- 
tion has been offered in Ref. [Sg- It was suggested that 
localised levels on the nanotube tip play the dominant 
role during the light emission. Therefore, measuring the 
noise power spectrum along with the luminescence spec- 
trum during FE experiments and comparing them with 
the prediction H5U|) would allow one to check the hypoth- 
esis put forward in Ref. j38j |. However, we would like 
to postpone the detailed discussion of this very extensive 
issue to a later publication and rather concentrate here 
on the ratchet effects. (We shall still give some general 
formulae for the noise power spectrum.) 



B. The non-interacting case 



We first discuss the case of non-interacting leads (and 
Ac = 0). Then the lead electron fields can be inte- 
grated out exactly. As a result we are then left with 
zero-dimensional problem of two fermionic levels. The 
corresponding Kcldysh GFs at zero-order in tunnelling 
are given by the following matrices: 



1 



{lU - Ai^2 



uj - Ai,2 - iVL^R sgn(w T V/2) i2TL,RQi~^ ± V/2) 

-i2rL,fle(t^ T v/2) -u) -f Ai,2 - iTl,r sgn{uj T V/2) 



(51) 



r 



+ 



FIG. 6: Diagrammatic representation of the Dyson equation 
1521 1. Solid lines represent the GFs of the level Ai and the 
dashed lines those of A2. Thick lines are the exact GFs. 

The energies Ai_2 are already renormalised by the gauge 
transformation and T r^l = Pr.liIr l where pr.l are the 
density of states (DOS) in the corresponding lead. In 
the absence of an external bias voltage V = ^l. The 
tunnelling term H\nr is quadratic so that we can write 
down the exact Dyson equation for the full GFs of the 
dot levels. In terms of the Keldysh contour ordered GFs 
it can be written as 

Di,2{t-t') = Dfl{t-t') + W'^ I dti dt2 D["l{t ~ ti) 

Jc 

(52) 



X Di'^lih - t2)Di,2{t2 - t') . 



The corresponding diagrams are depicted in Fig|Sl 
Disentangling the indices one obtains the following set 
of equations, 

DW = D^S'i^) - E mKl^2i^)DTSico) , (53) 

m— ± 

with the kernels 

kUu^) = -VF^ E mD^T^{u;)D^^l"''{u;). 

m— ± 



Equation system H53|l is linear and can easily be solved 
for all GFs. We need only the following ones (we drop 
the trivial energy argument uj): 



D 



1.2 



M"r'^it2-+[i-^iT]M°r 

detKi 2 
detKi 2 



(54) 



where 



deti^i,2 = (i-i^r2")(i + ^i+2^) + ^r:2^^it2" 



ii + if^p 



Evaluation of the current with help of (|44|l and (|54f) yields 
at zero temperature the following result: 



rV/2 

I = Gq dujT{uj). 

J -v/2 



(55) 



where the transmission coefficient is given by 



^^'^^ " [{uj - Ai)2 + al][{uj - A2)2 + al] - 2[{lo - Ai){uj - A2) - aRUL] + 1 ' ^^^^ 



From now on all energy variables are normalised to W pected that the optical coupling is much smaller than 
and dimensionless, a_R,L = ^r,l/W. However, it is ex- 
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FIG. 7: Light induced current through a quantum ratchet 
system without the applied voltage as a function of the 
light frequency il at zero temperature. The parameters are: 
Ai,2 = and = 0.1. The lower curves correspond to 
a non-interacting system while the upper curves represent a 
dot coupled to LLs with g — 1/2. 

that to the contacting electrodes, so that we can expand 
(|56|l for large aji,L- As a result we obtain the transmis- 
sion coefficient in form of two superimposed Lorentzians 
at energies Ai 2 and widths Tb.^. This implies for the 
I-V, 



I{V) 



V/2 



(57) 



duj 



2[(c.-Ai)2 + ri][(c.-A2)2 + r|] 



The latter integral can easily be calculated in a closed 
form. The results are given the Appendix IVII In the 
case of symmetric electrode couplings, — Tl = F, the 
current induced by the photon absorption (we shall call 
it 'ratchet' current) decays according to ~ W'^T/ft^ at 



high frequencies. In the opposite infrared limit it varies 
linearly: 

w^n ^ A,- 



Ai - A2 A 



i=l,2 



r2 



In the intermediate regime the ratchet current has its 
maximal value around ~ Ai — A2, see Figd A slight 
shift towards higher frequencies is a result of the mutual 
level hybridisation due to tunnelling similar to that oc- 
curring in the case of a double dot. It turns out that 
coupling asymmetry does not result in any qualitative 
change in the induced current. The qualitative picture 
remains the same for all values of the optical coupling 
W. Depending on the sign of the external bias voltage 
the ratchet current enhances or suppresses the transport. 
This effect might be of immediate experimental relevance 
since it indicates the way the levels of the quantum dot 
(or orbitals in the case of molecular dot) are arranged 
with respect to contacting electrodes. 

As already discussed in Section ^ the knowledge of 
T{uj) enables not only to access the transport but also 
the zero-frequency noise properties using formula H29|l . 
For the calculation of the full frequency dependent noise 
spectrum the mere knowledge of T{uj) is insufficient be- 
cause one needs the transmission amplitudes [s^ . In such 
a situation one can use the relation H50() . Generally the 
optical coupling is expected to be relatively weak with 
respect to the lead-dot coupling so that we can calculate 
the emission spectra at the leading order using Eas. (|51|l : 



TT 



V/2 



(58) 



2 1-1 



X / &[(6-Ai)^+Fi]-^[(6 + c.-A2)^+r^] 

This integral can be performed analytically but the result 
is lengthy. So we shall discuss only the special case of 
high bias voltage, corresponding to the FE via localised 
states. The evaluation of (|58|) then yields 



(Ai + A2 + c^)2 + (rL-rfl)2 



TT [(w-KAi-A2)2 



ri 



^l? - 2Firl, 



r 



Since the emission process lowers the energy of the 
radiating system the actual spectrum is Q{—U!)S^ {lo). 
Then this result describes a Lorentz-shaped peak around 
A2 — Ai , which is at odds with the experimental finding 
of Ref. J where a superposition of two Gaussian peaks 
has been detected. There are three possible reasons for 
this discrepancy, we list them in the order of relevance. 
First of all, the measuring apparatus can superimpose 



its own sensitivity curve, which is usually of a Gaussian 
shape, over the actual luminescence spectrum. Secondly, 
since the localised states on the nanotube tip could ex- 
ist on the ends of free dangling C ~ H bonds the finite 
temperature can lead to oscillations of the energy levels 
Ai^2, which then could become normally distributed. An- 
other possible reason could be the electronic correlations 
inside the nanotube. However, there is no a priori argu- 
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ment why their influence can result in Gaussian-shaped 
spectra. 



C. Interacting LL leads at g = 1/2 

Now we turn to the case of interacting leads. The 
Hamiltonian still contains the terms H^, He, the new 
tunnelling contribution (|41|) and the photon coupling 
term Hy/ after the gauge transformation (|42|l . As usual, 
we apply the transformation H32|l and in order to access 
the Toulouse point we again set the Coulomb coupling to 
Ac = Stt. In the language of the new fermions defined by 
(O we obtain then the following Hamiltonian (Ai_2 are 
assumed to be shifted by ±57/2): 



H = Ho[ij±] + ^ A,4d, 



1=1,2 



+ Ul 4'/'- (0) + IR 4^1 (0) + Wd\d2 + h.c. 



We proceed in the spirit of Section IIIII and derive the 
equations of motion for the participating operators, 

{idt - Ai)di^Wd2+"fLi'^iO), 

(idt " A2)d2^Wdi+jR'4'U0), (59) 
i[i/.U0+) - ^_(0-)] -7idi+7flVl(0), 

where 0^ is again positive (negative) infinitesimal. Act- 
ing with |yp [see Ea. ()36|) ] on the both sides of the third 
of Eqs. (|59|) we then use the first two ones in order to elim- 
inate the dot operators. As a result of this procedure we 
obtain an equation containing only operators, 

i|rnv^l(o+) - V-(0")] (60) 

^ tY*[W^R^pU0) + lLitdt - A2)V'-(0)] 

+ jRY[W^L^UO)+lRHdt - A2)^_(0)] . 

At this stage we again can make use of the decomposition 
(|38|l and reduce the relation (|60|) to exactly the form 
given in (|39|l with the same E [see definition (|40|) ] but 
different 7 and f3±, 

7 = Wy/^^]^{uj^ + AiA2-W^) , 
P+ - {aR[{uj + Ai){oj + A2) ~ W^]{uj ~ Ai) 

- aLiiu; - Ai)iuj ~ A2) ~ W^]{u; + A2)] /2 , 
P^iuj) = /3+(-^), 

where we again normalised all energy variables to F and 
defined aR,L = 7^ l/E. With these conventions the total 
dot transmission coefficient is still given by the formula 
© and the corresponding I — V hy Eq.©. 

In the case of a weak dot-leads coupling and the sym- 
metric level configuration (Ai^2 — ±const), the ratchet 
current does not show any significant change in compar- 
ison the the non-interacting case, see Fig[7| The only 



difference is a slightly higher and wider maximum. How- 
ever, as soon as the leads couple to the dot stronger than 
the dot levels among themselves an additional local min- 
imum emerges at low energies, see Fig|Sl This effect oc- 
curs only in the interacting system. More interesting 
features arise in a system with weak optical coupling, see 
Figini For positive fl (which in our picture corresponds 
to a situation favouring the ratchet effect) the distin- 
guished peak at Ai — A2 which exists at F <C splits in 
two (which is actually a more relevant parameter range 
from the experimental point of view). This effect does 
not occur in the non-interacting situation. The origin of 
the total four peaks can be traced back to the presence 
of the LL zero-bias anomaly in vicinity of the Fermi en- 
ergy. Concentrating only on the left half of the system - 
the left electrode with the level Ai, one can show in the 
lowest order in tunnelling that the correction to the level 
spectral function (which at 7^ = is a delta function) is 
given by 



6A{lo) 



{io-A^Y 



where v{u:) ^ [w/tJcl^^^^^ is the DOS of a half-open LL 
with a bandwidth lUc- At (7 = 1/2 SA{u!) possesses two 
maxima at ±Ai. This basic structure persists in all or- 
ders of 7l with some corrections to the maxima positions. 
Exactly the same thing happens in the right part of the 
system thus accounting for the total of four peaks seen 
in the ratchet current. 



V. SUMMARY 

To conclude, we presented exact solutions of the single- 
and two-state non-interacting QDs coupled to interacting 
(and non- interacting) electrodes. In both situations the 
corresponding Hamilton operators can be brought to a 
quadratic form, and thus solved exactly, at the special 
interaction strength g = 1/2 (the Toulouse point). 




FIG. 8: The same plot as in Fig|7|for the symmetric inter- 
acting system for different lead-dot coupling strengths. 
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FIG. 9: Ratchet current for a symmetric dot for different 
values of lead coupling. All energies are measured in units of 
|Ai|, Ai = -1, Aa = 4 and = 0.1. 

In the first setup, the single-state QD, we recovered all 
results previously obtained via scattering formalism with 
the help of the non-equilibrium GF approach. We have 
shown, that exactly at the resonance such model is fully 
equivalent to a simple tunnelling problem between two 
interacting leads at a reciprocal interaction parameter 
g = 2. Furthermore, we derived equations for all possible 
GFs and applied them to calculate the noise properties of 
the system which cannot be accessed by means of scat- 
tering formalism. Concentrating on the zero frequency 
noise properties we discussed the details of the voltage 
behaviour of the universal Fano factor vy- It turns out 
to interpolate between the non-interacting value at high 
voltages and unity in the asymmetric case (zero in the 
symmetric case) in the limit of small bias. In the asym- 
metric case possesses a minimum at some V* , which 
is absent in non-interacting systems. 

In the case of many levels on the dot the question of 
the conductance properties in the valley between the res- 
onant tunnelling peaks is important. In order to shed 
light on this, we considered a model of a two- level QD, 
in the case when each of them is coupled to both leads. 
This additional feature does not destroy the solubility of 
the system at the Toulouse point and we succeeded in 
applying the scattering formalism in this situation. As 
expected, the system turns out to show resonant trans- 
port signatures as long as the couplings are symmetric 
and one of the levels is tuned to the Fermi energies in 
the electrodes. If the gate voltage tunes the system into 
the valley between two peaks, G{T) still keeps its on- 
resonance temperature dependence with a renormalised 
pre-factor. However, turning on tunnelling between the 
two dot levels changes temperature dependence of the 
resonant transmission peaks in a non-trivial way. 



In the third model each level is assumed to be cou- 
pled only to one of the electrodes and the transport is 
supposed to be accompanied by absorption or emission 
of photons. It turns out that the current through the 
system can flow even in the absence of any the bias volt- 
age (ratchet current). In the non-interacting case the full 
non-linear I — V as well as the absorption and emission 
spectra can be easily calculated via the Keldysh diagram 
approach. We concentrated on the so called ratchet cur- 
rent induced solely by the electromagnetic irradiation in 
absence of any voltage sources. As expected, the depen- 
dence of the ratchet current on the light frequency pos- 
sesses a clear maximum at the energy difference of the 
dot levels. As soon as we take correlations into account, 
the picture changes considerably as the peak splits in two 
with a pronounced minimum between them. The origin 
of this suppression can be traced back to the zero-bias 
anomaly in the DOS of the interacting systems. Such 
spectacular effects make the ratchet current measure- 
ments an invaluable instrument for studying interacting 
QD structures. In addition we derived an exact relation 
between the absorption (emission) spectra and frequency 
dependent noise power spectrum. It has important impli- 
cations in the luminescence accompanying field emission, 
which is believed to occur during cold electron emission 
from carbon nanotubes. 

The key quantity, which generates an energy scale at 
which most of the predicted effects take place, is the 
lead-dot coupling F. In the typical experiments made 
on semiconducting QDs F ranges between 0.1 — 1 /leV, 
which corresponds to temperatures around 1 — 10 mK 
[40| . In the most current experiments conducted on con- 
tacted molecules the coupling strength is expected to be 
even smaller HQ. For the results of SectionslmlandlnBl 
to be accessible in the experiments it is not necessary to 
go significantly below the temperatures T ^ F. On the 
contrary, the experimental observability of all other pre- 
dicted phenomena, and especially of the ratchet effects, 
depends crucially on the ability to either lower the tem- 
peratures beyond the T ^ F mark or build devices with 
high enough F ,45j . 
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VI. APPENDIX A 

The result of integration in (|57|) is 
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